Electric field effects in STM imaging 
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We present a high voltage extension of the Tersoff-Hamann theory of STM images, which includes 
the effect of the electric field between the tip and the sample. The theoretical model is based on first 
principles electronic structure calculations and has no adjustable parameters. We use the method to 
calculate theoretical STM images of the monohydrate Si(100)-H(2xl) surface with missing hydrogen 
defects at —2 V and find an enhanced corrugation due to the electric field, in good agreement with 
experimental images. 



I. INTRODUCTION 

First principles electronic structure calculations have 
become an important tool in interpreting STM experi- 
ments. Calculations of theoretical STM images are often 
based on the theory by Tersoff and Hamann Q, which 
states that the the STM current is proportional to the 
local density of states(LDOS) of the sample. In this the- 
ory it is assumed that the potential is flat between the 
tip and the sample, and the vacuum level given by the 
sample work function. However, for relatively high biases 
(> 2 V), which are often used in STM experiments on 
semiconductor surfaces, the electric field strength in the 
tunnel region can be relatively high and must be included 
in the theoretical model ||. 

In this paper we extend the Tersoff-Hamann formalism 
to include the electric field in the tunnel region, and apply 
the theory to calculate the corrugation of a single missing 
hydrogen defect on the monohydride Si(100)-H(2x 1) sur- 
face. We find that the corrugation is strongly increased 
by the electric field, mainly due to polarization effects 
and partly due to changes in the tunnel barrier. 

The organization of the paper is the following: In sec- 
tion II we present the basic theory for calculating field 
dependent STM images, and in section III we show how 
the electric field effect can be included in the first prin- 
ciples calculation. In section IV we apply the formalism 
to calculate the corrugation of a missing hydrogen de- 
fect on the monohydrate Si(100)-H(2x 1) surface and in 
section V we conclude. 



II. THEORY 

In this section we present the basic theory for calcu- 
lating field dependent STM images. The derivation will 
follow Chen || closely. Figure 1 shows the tunnel junc- 
tion between the tip and sample. Using the modified 
Bardeen approach |3|, the tunnel current is given by 

1 = - eVh) - fi^u 2 ^, - a (i) 



where the integral in Eq. (2) is over any separating sur- 
face £ lying entirely within the vacuum region separat- 
ing the two sides. The sample bias, Vb defines the differ- 
ence between tip and sample Fermi levels, and /(e) is the 
Fermi function. The modified samplc(tip) wave functions 
(Xv) are solutions to the Schrodinger equation with 
modified sample(tip) potential U s (Ut). These potentials 
are given by the tunnel potential U upto the separating 
surface, S, and are equal to the vacuum level beyond 
the separating surface, thus U = U s + U t and U s U t = 0. 
The gradient of U in the tunnel region determines the tip 
induced electric field, E = VU/e. 
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FIG. 1. The STM tunnel region and tunnel potential. 

We assume that the tip is grounded and the vacuum 
level determined by the tip workfunction. Furthermore, 
we assume that the tunnel current is due to a single atom 



f 



at the tip apex. We place the separating surface, S, just 
outside the atomic radius, R, of this atom (see Fig. 1). 
The position of the tip we denote by r , and the tip sam- 
ple distance is given by d = r<j — R. Since the tip potential 
equals the vacuum level beyond the separating surface, it 
is straightforward to expand the modified tip wave func- 
tions in real spherical harmonics, l^ m , and obtain 



X 



„(r) =Y J C? m h{K t r)/k l {K e R)Y l m {r), 



CL = xARr)Yi m (r)dr, 



in 



k, = ^2m(0* + eVb + e£-e)/fi, 



(3) 

(4) 
(5) 



where r is the distance from the tip atom, ki the 
spherical modified Bessel functions, K e the inverse de- 
cay length of the electron states in vacuum, and e s F the 
Fermi level of the sample. Following Ref || we observe 
that the tip Green's function, defined by (— fi. 2 V 2 /2m + 
U t )G(r) = Aw6(r), is related to the tip wave-function by 
/c/(K £ r)yj m (f) = aiAi m G{r), where the differential 

operators Ai rn are defined i n Tab le 1, and t he coeffi cients 
ai defined by a.{ S:P:d } = 
now obtain the current 
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K l e +1 h(K e R) 



xD lm (e - eV b )de, 



(6) 



where we have neglected coherence between partial tip 
states and Di m (e) — J2u \^hn\ 2 ^( e ~ e ») are partial tip 
density of states per unit volume. 

We have calculated Di m for a single W atom on a 
W(110) surface using R = 3 bohr. We find that it is 
nearly independent of m, and average I dependent val- 
ues are D {s ^ d] (e) « {0.002e [ _ 10i20 ](e), 0.002e [0i20] (e), 
0.002Or_4 5 3 .51(e)} eV _1 bohr~ 3 , where the step function, 
®[a,b](%)i is one f° r x m the interval a < x < b and zero 
otherwise. Using these values we obtain the current 



,+eV b 



Y,Bi m (e-eV h )pi m (d + R,e)de, (7) 



metallic sphere of radius R, consistent with the spher- 
ical potential- well model of the tip used in Eq. (3) |Q. 

The main result of this paper, Eq. (7) , is a high voltage 
generalization of the Tersoff-Hamann expression jjj for 
the STM current. The main differences between Eq. (7) 
and the expression by Tersoff and Hamann are the inte- 
gration over the electronic states and the calculation of 
the sample wave functions in an external electric field. 
We also include higher angular tip states || , whereas the 
Tersoff-Hamann formulation is for an s-type state only. 
For the systems we have investigated we find that m > 
terms are more than one order of magnitude smaller than 
to = terms and can therefore be neglected. Of the 
to = states, we find that the / = state gives a con- 
tribution which is twice that of / > states. In the fol- 
lowing we will only consider the I = contribution, since 
we have found that this contribution best describes the 
experimental images we consider. However, we note that 
occasionally we see a change in the image contrast, which 
might be due to dominance of I > states for special tip 
geometries. 



/ m 



Bi m (e) 



A, 



0.007^e [ _ 10 ,20](e) 
0.02i? 4 (l + KeJ R)- 2 e [0 ,2o](X> 
0.02i? 4 (l + K£J R)- 2 e [0 ,2o](e) 
0.02i? 4 (l + K £J R)- 2 e [0 ,2o](e) 
0.03i? 6 (3 + 3k £j R + K 2 J R 2 )" 2 e [ _ 4 .5,3.: 



0.03i? 6 (3 + 3k c R + K 2 i? 2 )~ 2 e 



[-4.5,3.5] 



0.037? 6 (3 + 3k,R + K 2 i? 2 )" 2 e [ _ 4 .5,3.5] 
0.03.R 6 (3 + 3k,R + K 2 i? 2 )" 2 e[_ 4 .5,3.5] 
0.03.R 6 (3 + 3k,R + k 2 j R 2 )- 2 9 [ _4.5,3.5] 
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TABLE I. Definition of the parameters Bi m and differential 
operators Ai m used in Eq. (7). The step function, ®i a u(x), 
is one for x in the interval a < x < b and zero otherwise. 



where distances are in bohr, energies in eV and current 
in Amperes. Parameters Bi m are defined in Table 1. 
The main quantity is the sample LDOS, pi m {~x.,e) = 

53u I^WVv( x )l 2< ^( e — e ju)- The wave functions V-V are 
calculated in the external electric field from the tip, and 
we approximate this field by a planar electric field of 
strength E. For a given tip-sample distance d, the field 
strength is determined from the equation 



U s (d,E) = <f> t + eV h + , 



F ■ 



(8) 



where U s is the effective sample potential in planar field 
E. This equation assumes that the tip behaves as a 
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III. CALCULATION OF SURFACE ELECTRONIC 
STRUCTURE IN FIELD E 



calculations are based on 
using the Generalized 



The electronic structure 
density functional theory 
Gradient Approximation of Ref. |6j] for the exchange- 
correlation energy. Ultra-soft pseudo potentials in- 
cluding the nonlinear core correction j8| are used to de- 
scribe hydrogen and silicon and the wave functions are 
represented in a plane-wave basis set with kinetic en- 
ergy cutoff 20 Ry. At distances larger than 4 A from 
the surface, the wave functions are obtained by outward 
integration using the average effective potential perpen- 
dicular to the slab ||. 

In the following we will consider the monohydrate 
Si(100)-H(2xl) surface, which we model by a (2x1) slab 
with 12 layers of silicon atoms, and a vacuum region of 
10 A. We apply an external electric field to the surface, 
by inserting a dipole layer in the middle of the vacuum 
region ju^]. The effect of mobile carriers is introduced 
by fixing the atoms on the the back surface of the slab 
in their bulk positions ]li]j . This gives rise to half-filled 
surface states in the middle of the band gap, and this 
surface is therefore metallic. Depending on the direction 
of the external field the surface states accept or donate 
electrons. To obtain the field dependence of the wave 
functions we calculate the wave functions for two fields, 
Ei, E2 , which bound the field range in the experiment. 
The wave functions at a given field, Ei < E < E2, are 
then obtained by logarithmic interpolation between the 
wave functions at Ei and E2. 

Figure 2 shows the effective one-electron potentials for 
calculations with external fields of E = 0.6 V/A and 
E = —0.8 V/A. While the potential on the back surface 
is the same for both fields, the potential on the mono- 
hydrate(front) surface bends upwards for positive fields 
and downwards for negative fields, respectively. 

The inset in Fig. 2 shows the band bending in the slab 
calculation compared with the band bending in a n-type 
sample with Nd — 10 18 cm" 3 . The band bending cal- 
culation is based on standard band bending theory with 
non-degenerate statistics. |l^] For the n-type sample the 
Fermi level is close to the conduction band, and we have 
the well known depletion and inversion for positive fields. 
In the case of the slab calculation, mobile carriers are 
simulated by the half-filled dangling-bond states on the 
back surface, and since these states fix the Fermi level in 
the middle of the band gap, the band bending is nearly 
symmetric in the field. The experimental Fermi level is 
given by O 



Nn 

e s F =e q /2-0A9kT + kTlog(—) 

m 

4.73 x 10" 4 T 2 

£g = L17 ~ r + 636 eV ' 

n,i = 10 16 r5 e e f/ 2fcT cm" 3 , 



(9) 
(10) 
(11) 



where e g is the bandgap, rii the number of intrinsic car- 
riers and T the surface temperature. 

To correct for the difference in Fermi-level and band 
bending, $, between the experiment and the slab model 
we shift the STM voltage in the slab model by AV, + = 
$ + e g — €p — (<fr + e g — e s F ) at positive bias, and AV b ~ = 
<f> — e s F — ($ — e s F ) at negative bias (values with tilde are 
slab quantities). In this way we obtain that the energy 
window of electronic states which contributes to the cur- 
rent is the same in the slab model as in the experiment. 
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FIG. 2. The solid lines show the average effective potential, 
v e g, along the z axis(perpendicular to the slab) for external 
fields of 0.6 V/A and -0.80 V/A. The zero of the z axis is taken 
at the position of the first layer Si atoms, and the Fermi level 
defines the zero effective potential. The inset shows the band 
bending in the slab calculation(dashed line) compared with 
the band bending of a n-type sample with No = 10 18 cm -3 
at room temperature(solid line) lid] . 



IV. THE STM CORRUGATION OF 
SI(100)-H(2xl) 

In Fig. 3 we show a typical STM filled state image 
of the monohydrate Si(100)-H(2x 1) surface. The bright 
vertical stripes originate from the hydrogen passivated 
silicon dimer rows, and the white spot originates from 
a silicon dangling bond due to a single missing hydro- 
gen defect. Below the image we show the corrugation 
across the defect (solid line), compared with the simulated 
STM image of an s-state tip including field effects (dashed 
line) and without field effects (dotted line). In the range 
-5 A< x <10 A the theoretical curves were obtained us- 
ing a c(4x4) cell with a single missing hydrogen defect, 
and outside this range using a (2x1) cell. We see that 
the corrugation of the defect is well described in the field 
dependent calculation, while it is less than half the value 
when the field is not included. The larger corrugation 
of the defect in the field dependent calculation is mainly 
due to polarization of the dangling bond. Away from the 
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defect the corrugation in both calculations is less than in 
the experiment. Calculations using p- or (i-state tips do 
not give better agreement, and we suggest that the larger 
experimental corrugation away from the defect might be 
due to thermal vibrations on the surface. 
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FIG. 3. Filled state STM image of a single hydrogen defect 
on the Si(100)-H(2xl) surface, recorded with V = -1.6 V 
and 1 = 1 nA. The plot shows the corrugation across the 
defect (solid line), together with the theoretical corrugation 
including electric field effects (dashed line) and without field 
effccts(dotted line). 



V. CONCLUSIONS 

We have presented a high voltage extension of the 
Tersoff-Hamann model of STM images, which includes 
the electric field between tip and sample. We have ap- 
plied the model to describe the corrugation of a single 
missing hydrogen defect, and find good agreement with 
experiment when field effects are included. At low volt- 
ages Vb < 3 V, the field induced change of the corru- 
gation is mainly due to polarization. For higher voltages 
Vb | > 3 V, as used in many atom manipulation experi- 
ments, the field also has a pronounced effect on the tun- 
nel barrier, and we hope that the present work may prove 
useful for the analysis of such experiments. 



[1] J. Tersoff and D. R. Hamann, Phys. Rev. B 31, 805 
(1985). 

[2] H. Ness, A. J. Fisher, and G. A. D. Briggs, Surf. Sci. 
380, L479 (1997). 

[3] C. J. Chen, Introduction to Scanning Tunneling Mi- 
croscopy (Oxford University Press, New York, 1993). 

[4] P. Hohenberg and W. Kohn, Phys. Rev. 136, B864 
(1964). 

[5] W. Kohn and L. J. Sham, Phys. Rev. 140, A1133 (1965). 

[6] J. P. Perdew et al, Phys. Rev. B 46, 6671 (1992). 

[7] D. Vanderbilt, Phys. Rev. B 41, 7892 (1990). 

[8] S. G. Louie, S. Froyen, and M. L. Cohen, Phys. Rev. B 

44, 8503 (1991). 
[9] J. Tersoff, Phys. Rev. B 40, 11990 (1989). 
[10] J. Neugebauer and M. Scheffler, Phys. Rev. B 46, 16067 
(1992). 

[11] P. Kratzer, B. Hammer, F. Grey, and J. K. N0rskov, Surf. 

Rev. Lett. 3, 1227 (1996). 
[12] R. Seiwatz and M. Green, J. Appl. Phys. 29, 1034 (1958). 
[13] S. M. Sze, Semiconductor Devices (Wiley, New York, 

1985). 



4 



